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ABSTRACT 


An experimental program to obtain uniform inlet flow to 
the test blading in a large cascade facility designed to use 
inlet turning vanes, and to measure the conventional blade 
element performance, is described. Attempts to reduce non- 
uniformities (+1% in velocity) using screens were unsuccess- 
ful and so abandoned. Preliminary DCA blade element 
performance data were obtained without screens at one in- 
cidence angle before aero-mechanical problems with the inlet 
guide vane assembly curtailed testing. The blade surface 
pressure distribution at the one test condition compared 
very favorably with the distribution predicted using the 
NASA computer code QSONIC. Recommendations were made that 


would avoid the aero-mechanical problems encountered. 
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LIST OF SYMBOLS 


AVDR Axial velocity-density ratio 

Cp Coefficient of pressure at the inlet 
Cp, Coefficient of pressure at the outlet 
CPopartec Coefficient of static pressure rise 
Cc Blade chord (inches) 

D Diffusion factor 

i Incidence angle (degrees) 

P Pressure (in. H,0) 

Q Dynamic pressure (in. H,0) 

+ 3 Temperature (PR) 

X Non dimensional velocity 


8 Air angle, measured in the cascade midspan plane 
with respect to the axial direction (degrees) 


Y Stagger angle 
oO Solidity (C/S) 
w Loss coefficient 


Subscripts 

amb Ambient 

P Pressure 

PLENUM Plenum (supply) 

Ss Static 

wl North wall, lower plane 
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I. INTRODUCTION 


The need for lightweight, fuel efficient gas turbines 
that are capable of developing large amounts of thrust or 
power has motivated a continuing drive to obtain more ac- 
curate predictions of the flow through turbomachinery. 
Cascade testing of blade rows has, in the past, been a logi- 
cal and relatively inexpensive way to learn more about the 
phenomena involved in the flow through compressor and tur- 
bine stages. It is required more today in order to verify 
two dimensional and near-two dimensional analysis codes for 
flow through cascades. Such testing also provides two- 
dimensional blade element performance data which, in the 
absence of reliable analytical predictions, are required in 
the design of compressors and turbine stages. Reference 1 
describes how cascade measurements are obtained using a cas- 
cade wind tunnel and then used in the design process. 

Before subsonic cascade wind tunnel data can be accepted 
as being valid, the flow conditions must meet three require- 
ments. These criteria are discussed in detail in Refs. l, 2 
and 3. First, any disturbance in the airflow should be 
caused by the test blades; that is, the inlet flow to the 
test section must be acceptably uniform. 

Secondly, the measured flow characteristics should, 


ideally, be independent of spanwise position along the test 
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blades. The flow, ideally, should be two dimensional. 

Duval [Ref. 3] demonstrated that excellent flow conditions 
could be achieved in the Naval Postgraduate School Turbo- 
propulsion Laboratory (NPS/TPL) Subsonic Cascade Wind Tunnel 
using test blades with an aspect ratio of approximately two. 
The absence of suction along the walls results however, in 
some degree of streamline contraction which is measured in 
terms of an Axial Velocity Density Ratio (AVDR). 

The third requirement which must be satisfied is the 
periodicity of the inlet flow to the test section and of the 
outlet flow. Within one chord length of the leading edges 
of the test blades an upstream perturbation occurs as the 
streamlines adjust to negotiate the blade passages. Since 
the rectilinear cascade is simulating an infinite cascade of 
blades, the flow characteristics should be the same at cor- 
responding axial and blade-to-blade positions within each 
blade passage. This same condition should be true at any 
measurement plane downstream of the test blading. 

As described by Rose and Guttormsen [Ref. 4] several 
unique features were incorporated into the design of the 
NPS/TPL Cascade to ensure a two-dimensional and periodic 
flow at the test blading. Initial evaluations of the fa- 
cility were conducted and reported in Refs. 3, 4, 5 and 6. 
Work by Moebius [Ref. 6] involved modifications to the tun- 


nel plenum chamber which established satisfactory uniform 
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flow at the exit of the bellmouth contraction into the test 
section. 

In order to maintain an aspect ratio close to 2.0 ata 
solidity of 1.67 Cina [Ref. 7], following the work of Duval 
{Ref. 3], used a cascade configuration of 20 blades with 3 
inch spacing. Cina conducted a program of tests of DCA 
blading at five (5) different air incidence angles. With 
this cascade configuration, Cina found that the inlet flow 
to the test section was uniform in direction and of uniform 
static pressure, but with an imposed variation in velocity 
and stagnation pressure resulting from the wakes of inlet 
guide vanes. Although excellent periodicity was found over 
pairs of test blades, departure from strictly periodic con- 
ditions were detected from one blade passage to another. 
Cina explained this condition as being the result of the 
inlet guide vane wakes being separated at two inch intervals 
and entering a test section configured with a three inch 
blade spacing. Because of these flow conditions, Cina con- 
sidered his results to be preliminary. 

As a result of these findings the Cascade Wind Tunnel 
was modified so that inlet guide vanes were provided at one 
inch intervals. The object of the study reported herein 
was to obtain blade performance data on Cina's cascade with 
good periodic and uniform flow conditions. A necessary 
condition was to obtain agreement in the results for blade 


forces evaluated from surface pressures and from a momentum 
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balance. A second objective was to compare measured blade 
surface Mach numbers with the results of code calculations. 

At the outset, it was first necessary to carry out an 
extensive testing program to verify the new inlet guide vane 
section and the effect the new spacing had on flow uni- 
formity and periodicity. It was found that the uniformity 
of dynamic pressure was improved with the inlet guide vanes 
spaced at one inch intervals. Attempts were made to further 
improve the flow by the use of (various) wire screens placed 
downstream of the inlet guide vanes. These methods proved 
unsuccessful for the range of parameters tested and in fact 
aggravated the situation. 

Cina's testing of the Double Circular Arc blading was 
repeated without screens and with the Cascade Wind Tunnel 
configured with the modified inlet guide vane arrangement. 
Limited measurements were obtained before aero-mechanical 
problems with the new IGV arrangement, at the higher tunnel 
speeds, were encountered. 

The overall purpose of the testing program initiated by 
Cina was to obtain data with which to verify design optimi- 
zation computer codes developed by NASA. Towards this goal 
a fast, reliable computer analysis code (QSONIC) for cal- 
culating the flow field about a cascade of arbitrary 2-D 
airfoils was obtained from NASA. The code was adapted and 
modified to run on the Naval Postgraduate School's IBM 


370/3033 computer. 
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The program QSONIC was developed by NASA to overcome the 
Mach number limitations of the earlier program TSONIC [Ref. 
8]. QSONIC is described in Ref. 9. Procedures for using 
the program QSONIC at the Naval Postgraduate School's com- 
puter facility are given in Appendix D. The procedures are 
documented for the case of the DCA blading in the NPS/TPL 
cascade wind tunnel. A program listing is included to docu- 
ment changes made to the code in order to adapt to the oper- 
ating system of the NPS computer. 

Preliminary results show that experimental measurements 


and code predictions are in very good agreement. 


i? 





II. FACILITY DESCRIPTION AND MEASUREMENT APPROACH 


A. SUBSONIC CASCADE WIND TUNNEL 

The Naval Postgraduate School's Rectilinear Cascade 
Facility is shown in Fig. 1. A description of the facility 
as it was originally configured is given in Ref. 4. The 
test facility is an open cycle wind tunnel, designed for 
the purpose of testing cascades of axial-flow turbomachinery 
compressor or turbine blades. The unique design of the test 
section ensures that the airflow paths from the inlet guide 
vanes to all of the blades of the cascade test section are 
of equal length. This particular design was intended to 
eliminate the problems found in other cascade wind tunnels 
caused by having wall boundary layers of different thick- 
nesses entering the cascade at different points. 

As a result of the work reported in Ref. 1, two fine 
mesh screens were installed at the bellmouth entrance to 
improve flow stability. A follow-on study into the cascade 
performance was conducted by Bartocci and is reported in 
Ref. 5. As a result of Bartocci's findings, plenum turning 
vanes were installed to direct plenum inlet air towards the 
bellmouth entrance and to decrease the total pressure fluc- 
tuations present at the bellmouth entrance. Figure 2 shows 
the configuration of the plenum chamber as modified by Bar- 


tocci. Reference 6 describes work by Moebius that resulted 
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in further modification to the plenum chamber in which the 
original contraction was changed to two two-dimensional con- 
tractions in series. After this modification, acceptably 
small variations in velocity and flow angle were measured 
at the inlet guide vane station. Figure 3 shows the in- 
ternal arrangement of the plenum chamber as modified by 
Moebius and as it was configured for the work presented 
here. 

Using the plenum configuration shown in Fig. 3, Duval 
[Ref. 3] found that the wakes from the inlet guide vanes 
were not mixed out at the lower measuring plane of the test 
section but gave a well defined periodic variation in the 
impact pressure. The peak-to-peak variation was +4% of 
dynamic pressure over two-inch periodic intervals. This 
condition was undesirable, but was tolerated while looking 
only to establish the values of parameters recauired to 
achieve two-dimensionality and periodicity. Since the 
inlet flow conditions were not uniform, mass averages were 
used to calculate properties at the inlet plane from probe 
measurements. 

In order to achieve a solidity of 1.67 anc aspect ratio 
of about 2, a blade spacing of 3 inches was required for the 
tests carried out by Cina [Ref. 7]. The tests showed unac- 
ceptable departures from blade-to-blade periodicity under 
conditions of high blade loading and the installation of 


additional guide vanes was recommended. The modification 
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to the inlet guide vane section of the tunnel resulting from 
Cina's findings is described in detail in Appendix A. 

In the present work, several tests were completed with 
the tunnel further modified by the introduction of wire 
screens between the inlet guide vanes and the lower plane of 
the test section. Appendix B describes the screen material 
and the criteria used to select the particular screens used 


in this study. 


B. INSTRUMENTATION 

The instrumentation used in the present study is that 
which is described in detail in Ref. 7. Twenty static pres- 
sure taps were located on the north and south side walls. 
The taps on the south wall were connected to a water mano- 
meter board so that the uniformity of the static pressure 
distribution of the inlet and outlet could be monitored 
visually. Additionally, one upstream tap on each wall and 
one downstream tap on each wall (near the centerline) were 
also connected to the Scanivalve so that these static pres- 
sures were recorded. 

Figure 4 shows the probe that was used for the upstream 
survey (at the lower plane). The probe was a United Sensor 
Corporation DA 125 probe, serial number A847-1, calibrated 
earlier at various Mach numbers and yaw angles in a cali- 
bration facility. The United Sensor Corporation DC-125-24- 


F-22-CD probe, serial number A981-2 (Fig. 5), which was used 
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at the upper plane was similarly calibrated. The charac- 
teristics of the probes were approximated analytically to 
facilitate automatic data reduction procedures. The cali- 
bration and application procedures were those given by Duval 
[Ref. 3]. Appendix B of Ref. 3 describes both the upstream 
and downstream probes in detail. The mounting and traversing 


mechanisms are described in Ref. 7. 


C. REFERENCE MEASUREMENTS 

Plenum chamber (supply) pressure and temperatures, and 
atmospheric pressure were recorded on each data scan. 
Plenum pressure was also displayed on a water manometer 
board. The total temperature in the test cascade was as- 


sumed to be the same as the plenum chamber temperature. 


D. TEST BLADING 

The double circular arc test blading modeled the midspan 
section of the stator of the compressor stage reported in 
Ref. 10. Coordinates describing the profile of the blading 
are listed in Table D-2. The leading edge and trailing edge 
are shown in detail in Fig. 6. A photograph of the center- 
most blade is shown in Fig. 7. 

The three blades centrally located in the cascade were 
constructed with surface pressure taps along the midspan 
section as shown in Fig. 8. The centermost blade had 19 
ports on each of the pressure and suction surfaces and one 


tap at the leading edge. The two blades adjacent to the 
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center blade had 3 surface pressure taps located on each of 
the pressure and suction surfaces. The surface pressure tap 


locations for the centermost blade are given in Fig. 9. 


E. DATA ACQUISITION, REDUCTION AND ANALYSIS 

Data were recorded, reduced, and plotted using the modi- 
fied Hewlett Packard HP-3052A Data Acquisition System shown 
in Fig. 10. Reference 11 describes the system in detail. 
The system incorporated a HP-9845A desktop computer as a 
controller, with all components connected on the HP-98034A 
HP-IB Interface Bus. A NPS/TPL HG-78K Scanivalve Controller 
with two 48 port Scanivalves allowed the programmed acqui- 
sition of probe and blade surface pressure measurements. 

The software used in the present study for acquisition, 
reduction and plotting of data were developed from software 
originally created by Duval and Cina. The programs are 
listed and described separately in Ref. 12. 

The uncertainties in the measurements are listed in 


Table l. 
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III. EXPERIMENTAL PROGRAM AND RESULTS 


A. PROGRAM OF TESTS 

The test program was in three phases. First, in order 
to verify the new inlet guide vane assembly, tests were con- 
ducted with no blading in the test section and with the 
upper and lower endwalls set parallel at 35° (design condi- 
tion), 30° and 50° with respect to axial. 

Secondly, tests were made of the effect of wire gauze 
screen materials in reducing non-uniformities in the flow 
entering the test section. Appendix B describes the type 
of screens used and how they were installed. 

The last phase of the test program was a continuation 
of the work initiated by Cina. Table II lists the cascade 
configuration tested. One test was completed successfully 
before aero-mechanical problems were encountered and testing 


was halted until the causes were analyzed. 


B. TEST PROCEDURES 
1. Cascade Adjustments 
In the first and second phases of testing, the same 
procedures were used to realign the cascade for each new 
configuration. The lower and upper end walls A aie set to 
the desired flow angle and the inlet guide vanes were set 


so that their trailing edges were approximately aligned with 
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the end walls. The flow was started, and the desired inlet 
dynamic pressure was set. All tests were run at an average 
dimensionless inlet velocity (X) of about .13, corresponding 
to an inlet flow dynamic pressure of 18 inches water. Be- 
fore recording data, the water manometer board was checked 
to ensure that the distributions of wall static pressures 

at the inlet plane and outlet plane were acceptably uniform. 
If required, the inlet guide vanes were adjusted to obtain 
uniform static pressure to within +0.5 inches of water. 

In the third phase of testing, initially the pro- 
cedures used by Cina [Ref. 7] were followed, namely: the 
lower end walls were set to the desired inlet air angle and 
the upper end walls were set approximately to the expected 
exit air angle. The inlet guide vanes were set very approxi- 
mately and the cascade was turned on and set to an inlet dy- 
NMamic pressure of 18 inches water. The upper end walls and 
the inlet guide vanes were adjusted in turn to obtain wall 
static pressure distributions upstream and downstream which 
were acceptably uniform. Using this procedure however it 
was found on occasion that the inlet air angle sensed by the 
probe at the lower plane at mid-span could be 2 or 3 degrees 
different from the setting of the end walls. 

The following procedures was subsequently adopted. 
The lower end walls were set to the desired inlet air angle. 
The upper end walls were adjusted to be "wide open", to form 


a diverging passage in which, when the cascade was turned 
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on (to an inlet dynamic pressure of 18 inches water), the 
flow was completely separated. The inlet guide vanes were 
adjusted to obtain the required inlet air angle on the chan- 
nel center line over the center 24 inches in the blade-to- 
blade direction. The upper end walls were then moved 
individually towards the vertical until the lower plane 
static pressure distribution was uniform and the upper plane 
static pressure distribution was acceptably uniform at a 
value close to atmospheric pressure. No readjustment of the 
inlet guide vanes was made. 
2. Measurements 

Probe surveys were carried out in the blade-to-blade 
direction at midspan at the upper and lower planes. [In the 
first and second phases of testing, data were taken over 
approximately 24 inches of the test section at intervals of 
0.25 inches. Also, in order to test the repeatability of 
measurements, repetitive samples were taken with the probe 
held fixed at midspan at the lower plane at the center, 10 
inches to the right and 10 inches to the left of center. 

During the third phase of testing, data were taken 


using the procedures established by Cina in Ref. 7. 


C. VERIFICATION OF INLET GUIDE VANE (IGV) ASSEMBLY 
The results of the first phase are presented (as shown 
in Table III) in Figs. 11 to 32. The results are arranged 


into groups. The first group (Figs. 11 to 14) are 
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measurements of tunnel conditions with the end walls set at 
35 degrees. Plots of conditions at the lower plane are fol- 
lowed by plots of conditions at the upper plane. 

Results for a wall angle of 30 degrees are given next 
(Fig. 15 to Fig. 18), followed by results for a wall angle 
of 50 degrees (Fig. 19 to Fig. 21). Data were taken over 24 
inches at the lower plane, and also at the upper plane at 
30°. At 50°, at the upper plane, only the center 12 inches 
were surveyed. 

The degree of repeatability of conditions in the wind 
tunnel from test to test (with no change in wall setting) is 
demonstrated by the results plotted in Figs. 22, 23 and 24. 

The last group of plots, Figs. 25 to 33, shows the de- 
gree of repeatability in the probe data from scan to scan. 
Data for these plots were obtained by holding the probe sta- 
tionary at midspan in three specific blade-to-blade loca- 
tions in turn and taking 50 repetitive scans of the channels 
normally recorded for survey profile data. The time inter- 


val for each scan was approximately 20 seconds. 


D. TESTING WITH WIRE GAUZE SCREENS 

The selection and installation of the wire gauze screens 
is described in Appendix B. 

The measurements obtained with the various screen con- 
figurations are given (as shown in Table IV) in Figs. 34 to 


45. The results are arranged in four groups. 
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The first group of plots (Fig. 34 to 40) give data ob- 
tained with the 16 mesh .0105 inch diameter wire screen in- 
stalled. Over a blade-to-blade distance from -1.0 to 7.0 
inches (Fig. 35) a peak-to-peak variation in velocity of 
about 1 percent was noted. This is slightly greater than 
the less than 1 percent (0.9 percent) variation noted over 
the same survey region without a screen installed (Fig. 12). 

The results shown plotted in Figs. 41 and 42 were ob- 
tained with two screens installed. One screen (16 mesh, 
-0105 inch diameter) was installed as discussed in Appendix 
B, while the second screen (2 mesh, .0400 inch diameter) was 
attached across the duct at the leading edges of the inlet 
guide vanes. 

The fourth group of plots (Figs. 43 to 46) show the re- 
sults for two different singhe screens. Probe survey data 
for these screens was taken only at the lower plane. The 
results shown plotted in Figs. 43 and 44 are data obtained 
with a 4 mesh, .041 inch diameter wire screen installed. 

The variation in velocity in the blade-to-blade direction 

was as much as +l.1 percent, neak-to-peak. The results 

shown in Figs. 45 and 46 are for a 5 mesh, .041 inch diameter 
wire screen. The variation in flow velocity was approxi- 


mately +1.5 percent, peak-to-peak. 
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E. PRELIMINARY TESTING OF DCA BLADES 

The results contained in Tables V to IX and Figs. 47 to 
60 are arranged in the following manner. 

The results shown plotted in Figs. 47 to 60 are divided 
into two separate groups. The first group (Figs. 47 to 57) 
contain results which exhibit the quality of the wind tunnel 
flow conditions. The second group (Figs. 58 to 60) shows 
the blade forces (and surface pressures) from survey data. 
In the first group of figures, results are presented first 
to examine the inlet flow uniformity (Figs. 47 and 48); 
second, to examine the outlet flow periodicity (Figs. 49 to 
53); and, finally, to examine outlet flow two dimensionality 
ieage. 5S¢ to’ 57}. 


All points are shown connected with straight lines. 
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IV. DISCUSSION OF EXPERIMENTAL RESULTS 


A. EFFECT OF INLET GUIDE VANE (IGV) MODIFICATION 

The probe survey data shown in Figs. 11 and 12 were 
taken at the lower plane, with the end walls set at 35 de- 
grees. A turning angle of 35 degrees corresponded to the 
"design point" of the inlet guide vanes, when the airflow 
from the plenum chamber was at zero angle of incidence to 
the leading edge of the IGV's. These two figures show that 
the inlet plane total pressure at midspan in the blade-to- 
blade direction had a peak-to-peak periodic variation of 
about +2 percent and therefore about a +l percent peak-to- 
peak variation in the velocity. Corresponding data from 
probe surveys at the upper plane (Figs. 13 and 14) show that 
periodic variations in total pressure were reduced to about 
25% of the value at the lower plane by the mixing of the in- 
let guide vane wakes. 

Tt can be seen in Figs. 15 to 21 that at "off-design" 
conditions for the IGV's (endwalls at 30° and 50°) there is 
a greater periodic variation in total pressure at the lower 
plane in the blade-to-blade direction than at the design 
point conditions. In Fig. 15 and Fig. 16, it can be seen 
that the periodic variations in total pressure are more pro- 
nounced with the flow from the plenum at a negative incidence 


angle to the IGV's (endwalls at 20° \. Except for the first 
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8 data points in Fig. 15, there is a well defined period of 
about 1 inch of travel. 

Figures 19 and 20 show the probe surveys conducted with 
endwalls at 50 degrees. At this off design condition the 
periodic variation in total pressure was considerably 
greater than the design point and the individual wakes from 
the inlet guide vanes were much less well defined. 

The repeatability of the survey was examined at wall 
angles of both 50° and 30°. Figs. 22 to 24 show that the 
non-uniformities in the flow conditions were reveated to 
(generally) better than 0.5 percent of total pressure. The 
question was then, to what accuracy could the individual data 
points be repeated in successive samples. This was examined 
at several probe positions and the results given in Figs. 25 


to 33 explain the departures in Figs. 22 to 24. 


B. EFFECT OF WIRE GAUZE SCREENS 

All testing with screens was conducted with the end walls 
and inlet guide vanes set to yield a flow angle of 35 de- 
grees. The data obtained with screens installed were there- 
fore compared with the data obtained without screens, shown 
in Figs. 11-14. The effect of the pressure drop across the 
screen on the pressure coefficient plotted in Figs. 34-35 
should be noted. With the first screen installed the drop 
in total pressure from plenum to the probe in the lower plane 


was about 10 inches of water (plenum pressure minus total 
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pressure measured by the probe at the lower plane). With- 
out the screen (at design conditions), the pressure drop 
from the plenum to the probe at the lower plane was approxi- 
mately 2.0 inches of water. Since Q-a¢ Was defined as the 
difference between plenum pressure and lower wall static 
pressure, the value of i with the first screen installed 
was about 28 inches water and without the screen installed, 
about 20 inches water. In comparing the peak-to-peak varia- 
tion in Pl seen in Fig. 11 with that obtained with the first 
screen installed in Fig. 34, the difference in the values of 


Q Must be considered. Examinations showed that the peak- 


ref 
to-peak variation in velocity remained at approximately 1% 
when the screen was installed. 

Figures 38 to 40 are plots of data obtained during a 
spanwise traverse of the probe at the lower plane. These 
figures show that the pressure drop through the screen and 
turning vanes was nearly uniform over approximately 8.0 
inches of the 10.0 inch span of the tunnel. 

In an attempt to generate upstream disturbances that 
Might trigger early boundary layer transition on the IGV's 
and increase the rate of mixing of the wakes, a second screen 
was attached to the leading edge of the IGV's. The results 
in Figs. 41 and 42 showed that this was not the case and in 
fact the second screen increased the magnitude of the non- 


uniformities at the lower survey plane. Measurements made 


with a single 4 mesh (Figs. 43 and 44) and a single 5 mesh 
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screen (Figs. 45 and 46) of similar blockage showed that 
neither screen influenced the flow in a particularly favora- 
ble manner. The 4 mesh screen caused the peak-to-peak 
variation in velocity to be about +1.1 percent, while the 
5 mesh screen caused the variation to be about +1.5 per- 
cent. This compared unfavorably with the variation ob- 
tained without any screen installed which was less than +1 
percent. 

It was therefore decided to proceed with measurements 


of the test cascade without using screens. 


C. PRELIMINARY TESTING OF DCA BLADES 
1. Inlet Uniformity 
The probe survey at the lower plane in Fig. 48 shows 
that the inlet plane total pressure at midspan varied in the 
blade to blade direction less than 0.5 inches of water, with 
no well-defined spatial period. This was an improvement in 
the inlet conditions found by Cina [Ref. 7: Fig. 16]. That 
the spatial period was not well defined agreed with the 
findings presented earlier in this report. The wall static 
pressure distribution (Fig. 47) showed small variations 
(less than 0.5 ins. water peak-to-peak at the lower plane, 
-4 ins. water peak-to-peak at the upper plane). 
2. Two-Dimensionality 
The data in Figs. 54 to 57 show that, at the down- 


stream plane, an area of (spanwise) nearly uniform conditions 
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existed near the centerline of the cascade. Reference 2 
points out that at higher loadings it is difficult to estab- 
lish a substantial spanwise area of uniform flow in the 
region near the suction side of the blade. This difficulty 
is evident in the data shown in Figs. 54 and 55 which show 
that only about 20% of the spanwise distance is acceptably 
uniform. It is noted that Cina also found reduced areas of 
uniform flow at this incidence angle; however, 30-40% of the 
spanwise distance was found to be acceptably uniform in his 
case. The difference could be the result of the reduced 
spacing of the IGV's and its effect on the side wall boundary 
layers. 

Figure 58 shows results for inlet and outlet flow 
angles and blade force vectors derived in two ways as shown 
in Appendix B of Ref. 7. These two methods are first the 
applications of momentum conservation to probe survey data 
and second, the integration of surface pressures measured 
over the blade area. Reference 2 points out that for truly 
two dimensional flow the blade forces derived from the two 
methods should be the same. As shown in Fig. 58 the magni- 
tudes and the directions of the two vectors representing the 
blade forces are in reasonable agreement. It is noted how- 
ever that at this particular incidence angle Cina [Ref. 7] 
Measured blade forces that were in total agreement in direc- 


tion but disagreed slightly in magnitude. The values of the 
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force magnitudes were about 1.5% lower than those measured 
by Cina. 
3. Periodicity 

As can be seen in Figs. 50 and 51 the total pressure 
and velocity qualitatively repeated fairly well over three 
central blade passages. Acceptably small quantitative dif- 
ferences are noted. There was also a small but measurable 
difference in the surface pressures on adjacent blades, as 
is evidenced in Fig. 49. 

4. Blade Performance 

Figures 59 and 60 are plots of the pressure and 
velocity distributions respectively over the centermost 
blade. These results compare favorably to those obtained 
by Cina for an incidence angle of ie eee 

Table IX contains the blade performance parameters 
deduced from the probe survey data listed together with the 
data obtained for corresponding test parameters in Ref. 7. 
While differences in two sets of data are evident, the dif- 
ferences are not large. It is noted that the value of the 
loss coefficient was only 10% lower than was measured by 
Cina, but the AVDR was less than 2% different from unity 
rather than the 6.5% measured by Cina. Further measurements 
need to be made, particularly in the light of the following 


discussion, before stronger conclusions can be drawn. 
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5. Aero-Mechanical Problems Encountered 

Fifteen cascade tests were made while evaluating the 
new inlet guide vane assembly and testing the wire screens. 
All runs were made without test blades installed, with a 
plenum total pressure of about 20 inches of water and with 
the upper and lower end walls parallel. No difficulties 
were encountered in establishing the desired flow conditions 
or in using the inlet guide vanes to arrive at a satisfactory 
distribution of wall static pressures. 

The first time the Cascade Wind Tunnel was set up 
with test blades installed to take data at an air inlet angle 
of 2.2", the tunnel operated normally and the test was com- 
pleted. (The data from this test were subsequently found to 
be highly suspect and are not reported here.) During the 
next test, with an inlet air angle of 42.4, the start-up 
appeared normal and previously established procedures were 
used to arrive at a satisfactorily uniform wall static pres- 
sure distribution. Tunnel operation appeared to be normal 
while taking data, but on shutdown a very noticeable high 
frequency vibration was encountered. Examination of the in- 
let guide vanes revealed that about 40% of the 60 blades 
were damaged. Damage included chips missing from the trail- 
ing edges, blades bent, cracks at the weld where the blade 
is joined to its support and indications that the suction 


side near the leading edge of one blade had been vibrating 
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against the pressure side near the leading edge of an ad- 
jacent blade. 

After the inlet guide vanes had been repaired and 
reinstalled extreme care was used at the beginning of the 
next test to adjust the inlet guide vanes, with two indi- 
viduals monitoring the movement of the adjustment mechanism. 
(IGV adjustment mechanism is described in detail in Appendix 
A.) A lack of stiffness in the mechanism was suspected as 
having been a contributing factor to the failure. 

One successful test was completed at an inlet flow 
angle of 42.4° and these data were discussed above. 

With the next cascade configuration set, at an air 
inlet angle of 45.9°, when the IGV's were adjusted after 
starting up, high frequency vibrations were again experi- 
enced. The wind tunnel was shut down and no further testing 
was attempted at plenum total pressures as high as 20 inches 
of water gauge. 

The difficulty encountered with the IGV's is not 
fully understood, however the lack of stiffness present in 
the actuation of the two separate rows of vanes is suspected 
of having allowed the problem to occur. Certainly, the pos- 
sibility of an aerodynamic flutter condition being present 
(due to the misadjustment perhaps) can not be ignored. It 
was noticed after the initial failure that the lead screw 
which adjusts the IGV's could be turned but the blades 


mounted from only one side would be caused to rotate. This 
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could lead to the trailing edge of one blade contacting the 
trailing edge of an adjacent blade and effectively closing 
the blade passages. 

Also, the holes in which the cylindrical shanks of 
the IGV's were held were found not to be uniformly machined. 
As much as 0.1 inches of movement at the tip of some vanes 
was possible while others could barely move. (The most 
seriously damaged blades were found in or adjacent to the 
larger holes.) 

The tendency for the mechanism to "hang-up" on one 
Side would be greater as the vanes became more highly loaded. 
It is noted that the IGV problem was encountered first when 
going to increased incidence angles with the compressor test 
cascade installed. In setting a constant plenum total pres- 
sure of 20 inches of water gauge, the static pressure in- 
crease across the test blades to a constant atmospheric 
pressure at the downstream side implies that a progressively 
increasing dynamic pressure was being generated out of the 
turning vanes. This can be seen in Table IX, where Q) for 


the test at By = 42.4 was 25 inches of water. 
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V. COMPUTATIONAL PROGRAM 


A. DESCRIPTION OF QSONIC 

The computational code, QSONIC, was developed by the 
staff at the NASA Lewis Research Center. This code is able 
to calculate the blade-to-blade flow conditions in turbo- 
machinery blade rows assuming inviscid flow but including 
streamtube convergence and radius change in the throughflow 
direction. QSONIC is flexible enough to allow the input of 
the appropriate boundary conditions to calculate the flow 
through the test blading in the Subsonic Cascade Wind Tun- 
nel. The program uses a fully conservative solution of the 
full potential equation combined with the finite volume 
method on a body-fitted mesh. QSONIC uses an artificial 
density imposed in the transonic region, if such a region 
exists, to ensure stability and the capture of shock waves. 

The analysis used by QSONIC is a combination of transonic 
analysis methods to calculate the flow conditions in the vi- 
cinity of a cascade of airfoils. A conservative form of the 
full potential equation is discretized at every point of a 
body fitted periodic mesh and a mass balance is calculated 


through the finite volume surrounding the point. The volume 





1the help and advice received from Charles Farrell at 
NASA Lewis R.C. in the process of adapting QSONIC to the 
NPS computer is gratefully acknowledged. 
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is corrected three dimensionally for any change in stream- 
tube thickness along a streamtube, if a quasi-3D solution is 
desired. Either elliptic or hyperbolic non-linear partial 
differential equations are used, depending on the local Mach 
number. 
The analysis used in developing QSONIC made the following 
assumptions: 
1) The airflow is inviscid and adiabatic. 
2) The airflow relative to the test blades is steady. 
3) Air is a perfect gas with constant specific heat. 
4) The airflow is isentropic and any discontinuities such 
as shocks are so weak that they may be approximated as 
isentropic jumps. 


5) There is no velocity component normal to the 
streamsurface. 


6) The airflow relative to a fixed reference frame (i.e. 
absolute velocity) is completely irrotational. 


Assumption 4 requires that the peak local relative Mach num- 
ber on a blade surface be 1.4 or less. The Mach numbers 
measured in test blades in the Subsonic Cascade Wind Tunnel 
would be well within this limit. However, this limitation 
would probably preclude the use of QSONIC for analysis of 
the flow field in the NPS transonic cascade wind tunnel. 
There are some combinations of blading geometry and flow 
conditions which cause unsatisfactory results to be gener- 
ated. For example, because of assumptions 1 and 6, sharp 
leading edges at high incidence angles (more that a few 


degrees) cause large velocity peaks in the blade surface as 
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the flow tries to turn from the stagnation point to the 
suction surface. 

Reference 9 gives a detailed description of QSONIC and 
the solution method used including the governing equations. 
Appendix D describes the operating procedures to use QSONIC 
on the Naval Postgraduate School's IBM 370/3033 computer. 
Appendix D also describes the input and output required as 


applicable to the Subsonic Cascade Wind Tunnel. 


B. APPLICATION TO THE TEST CASCADE 

Appendix D describes in detail the generation of the 
input required for QSONIC when applied to test blading in 
the Subsonic Cascade Wind Tunnel facility. In the present 
work one comparison of code calculations and measured data 
was made before testing was stopped. The comparison was 
for an inlet flow angle (85) oF 42.4". 

Tables D.1 and D.3 show the input cata generated. Table 
D.6 shows the flow solution output by QSOMNIC. The flow cal- 
culated on the blade surface, using a 15 by 97 mesh, was 
examined. Figure 61 is a plot of the calculated Mach number 
along the blade surface using two dimensional inputs. Fig- 
ure 62 is a plot of computed Mach number incorporating quasi- 
three dimensional effects. The method of incorporating 


quasi-three dimensional effects is explained in Appendix D. 
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C. COMPARISON OF CODE CALCULATIONS AND MEASURED DATA 

Table VII lists the data measured in the cascade wind 
tunnel. Cole Co2 and Xe] 2re defined in Appendix C. The 
surface Mach number distribution measured on the center blade 
is shown plotted in Fig. 63. 

For comparison, the computed two dimensional, comouted 
quasi-three dimensional, and the Mach number measured in the 
cascade wind tunnel are plotted together in Fig. 64. 

Excellent agreement between all three cases is seen. 

As would be expected, the greatest difference between mea- 


sured and calculated data is near the leading edge in the 


Suction side and at the trailing edge of the blade. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


Based on the results of the first part of the present 
Study, to evaluate the effects of the altered inlet cuide 
vane spacing on flow uniformity and periodicity, the follow- 


ing conclusions were drawn: 


1. With the inlet guide vanes operating at design, the 
peak-to-peak variation in velocity was about 11%, and 
there was a well defined spatial period of about 1 
inch. 


2. Operating the cascade wind tunnel in a configuration 
that requires the inlet guide vanes to be set to other 
than zero incidence resulted in peak-to-peak varia- 
tions in velocity greater than +1% and a spatial period 
that was less well defined. 


The second part of the study, to evaluate the use of 
wire gauze screens to further reduce the non-uniformities in 
the flow field, led to the following conclusions: 

.l. A 16 mesh screen with a blockage factor of .69 had a 
slightly aggravating effect in the variation in velo- 
city. The peak-to-peak variation in velocity with the 
screen installed was slightly greater than with no 
screen installed. This occurred at the expense of a 
pressure drop of about 10" of water across the screen. 

2. The use of screens with similar blockage but with 
larger mesh and larger diameter wires resulted in 
larger peak-to-peak variations in the velocity at the 
lower plane. 

The overall objective of the present study was to mea- 
sure the performance of the DCA test blading. Because of 


aero-mechanical problems encountered with the inlet turning 


vanes the performance of the blades was obtained at one 
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incidence angle only. The following was concluded from the 


limited test program: 


i. 


As a result of the reduced inlet guide vane spacing 
the variations in velocity and total pressure at the 
inlet plane were much less than those reported by Cina. 


Good periodicity was found from one blade passace to 
another. 


An acceptable region of spanwise uniformity (20-40% of 
blade span) was found at the downstream plane at the 
one test condition reported. However, this was less 
than was previously reported for the same incidence 
angle. 


The blade forces derived from the integration of sur- 
face pressure measurements and probe survey data were 
in close agreement in both magnitude and direction. 


The Mach number measured by surface pressure taps over 
the surface of the blade and Mach number calculated 
using the program QSONIC were in excellent agreement 
qualitatively and reasonable agreement quantitatively. 


The specific reasons for the aero-mechanical problems 
experienced with the inlet turning vanes have not been 
identified completely. 


Based on these conclusions and other observations, the fol- 


lowing recommendations are made: 


1. 


Use of the present inlet guide vane assembly and ad- 
justment mechanism for testing at inlet dynamic pres- 
sures higher than about 15 inches of water is unsafe. 
There are three possible solutions to this problem. 


a) Operate only within the dynamic pressure range of 
10-15 inches of water. 


b) Modify the new inlet guide vane assembly so that 
the vanes are supported at both ends on their axes 
of rotation. [Supporting the IGV's from both ends 
would prevent flapping vibrations of the (present- 
ly) cantilevered vanes. Such vibrations, when the 
vanes are supported alternately from opposite ends 
and the gaps are small compared to the chord, 
might lead to a potentially destructive flutter 
mode at particular flow velocities. ] 
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c) Replace the entire inlet turning vane section with 
one of entirely new design. 


The procedure should be adopted immediately of adjust- 
ing the vanes and walls of the cascade at lower values 
of the dynamic pressure before increasing the blower 
speed to the desired operating condition. 


More time needs to be spent to examine the flow field 
produced between the guide vanes and the test blades, 
and to establish the effects of the movement of the 
tail boards. The uniqueness of the flow field when the 
wall static pressure uniformity is used as a criterion 
of good inlet flow, needs to be examined by conducting 
repetitive tests at nominally similar test parameters. 
Only when the adjustment of the flow in the facility 
and the quality of the flow itself is fully understood 
should the measured blade performance data be accepted 
as final. 


The upper electrical yaw adjustment mechanism should 
be replaced with a manual system to greatly decrease 
the time required to achieve probe pressure (angle) 
balancing. 


Develop the computer code necessary to take advantage 
of the plotting data created by QSONIC. 


Modify the data acquisition and reduction software for 


the HP 9845 so that real time plots of blade perfor- 
mance parameters can be displayed. 
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TABLE I. MEASUREMENT UNCERTAINTY 


Item Description Method Uncertainty 
a gc nnn pense a ened incense 
x Blade-to-Blade dimension Position +01 ain, 
x = O in. West end Potentiometer 


x = 60 in. East end 








z Spanwise dimension Position =. 01, in, 
z= OQOin. North wall Potentiometer 
z= 10in. South wall on probe mount 

8, Inlet flow yaw angle Angle Potenti- +.2 deg. 


ometer on probe 
mount (hand 
adjustment) 





8, Outlet flow yaw angle Angle Potenti- t.5) deg. 
ometer on probe 
mount (motor 
driven adjustment) 








P a Plenum total pressure Static tap in Soke Ta. H,0 
P plenum chamber gauge 
V= 0 

P. Static pressure at the Calibrated pneu- eek in, 4,0 
test plane matic probe gauge 

Pn Static pressure at Static tap on =) - an: H,0 
=o 26.,- 9 = -16.23 in., North wall gauge 
z= 0 is. 

P m™ Atmospheric pressure Absolute Strain +3) Lis H,0 


Gauge Transducer 





P Pressure Scanivalve #204 in, H,0 
Transducer gauge 
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TABLE II. CASCADE CONFIGURATION FOR DCA BLADE TESTS 


Constant Parameters 


Number of Blades 20 

Spacing (Pitch) 3 inches 

Chord 5.01 inches 

Solidity DRE OW 

Thickness 7.0 percent of chord 
Camber Angle 45.72 degrees 
Stagger Angle 14.72 degrees 


Variable Parameters 


By 42.4 degrees 


i 5.3 degrees 
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TABLE III. SUMMARY OF MEASUREMENTS WITHOUT SCREENS 


Survey Plane Survey Direction Fig. Nos. Purpose 
Lower B=-B 1 & 2 
(24 inches) 


Upper B-B 3 <@ £4 
(24 inches) 


Lower B-B 15. & 16 
(24 inches) Flow Field 
Determination 
Upper B-B 17 & 18 
(24 inches) 
Lower B-B LO .& 20 
(24 inches) 
Upper B-B 21 
(12 inches) 
Lower B-B 22 & 23 
(24 inches) 
Survey 
Repeatability 
Lower B-B 24 
(24 inches) 
Lower Fixed Probe 25-27 
ino" o£. of ©) 
Point 
(on €) 28-30 Repeatability 
(10° RB. ot ¢.) 31-33 
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TABLE IV. SUMMARY OF MEASUREMENTS WITH SCREENS 


s fe) 
(By = 35} 
Survey 
Screen Survey Plane Direction 
16 mesh Lower B-B 
.0105 wire 
Upper B-B 
Upper Spanwise 
16 mesh Lower B-B 
.0105 wire 
+ 2 mesh 
ahead of IGV's 
4 mesh Lower B-B 
.041 wire 
5 mesh Lower B-B 
.041 wire 
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Fig. Nos. 


34 & 35 
3658 37 


38-40 


41 & 42 


43 & 44 


45 & 46 
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TABLE V. PROBE DATA, UPPER PLANE AT MIDSPAN (i = 5.3°) 


UHTAR FROM FILE WEREDe:Ti4 
BLADE TO BLADE TRAVERSE MID ZPAH 


UPPER PLANE 


Forme Loctin: Beta H- Oibar Feolbar Pr om ibar Mo Abe ae 
ie ai Ee ie ee ii ie Ee ie ee Dee Be ee De ee ee a ee Ee De ee Ba De ee ee te eae ee ee Bee eed 
1 =Fesi =n ae .3656 .2e14 B32 26475 
= -o.34 -.34 saose .2eG? 1841 6331 
3 =6. Sr =233 $276 256! 22336 Seiad) 
- -aeer =2556 . 3096 -2rar .1388 e145 
S -3.6f =-2.56 a te 22633 .as95 26425 
5 -5.245 25 2 23°36 20a . 9962 25958 
3 -3.23 tat ara s2vres 29916 Seheo tes: 
3 -5.84 “A579 5747 2626 Pee ea 25435 
3 -4.50 -1.78 ,seor ~2573 . 8352 .5436 
18 -4,.58 =1.73 .3f85 .2538 ~-898l 22935 
it -4.41 -1.79 - 3687 me oad 29955 2475 
l2 “4.15 =1.36 . 5658 senel 89394 2345 
13 -3. 96 =1.35 3629 . 2686 9987 26445 
14 =—Sute +i. 36 . 5462 aed 9 »1123 ,5a03 
15 -3.58 -1.36 .4553 2634 ies SBE 
is Sees eS =1..33 »4349 2618 5 BASS) Sali i 
17 =3. 89 -1.88 » 4239 vor e: 2455 ae ke ee 
13 -3.68 -3.16 ~4836 2°38 . 1668 Seis =) 
13 See $2.41 S467 2744 1168 Meh) i=, 
24 -2.548 =2. 493 sort? 2636 8902 6461 
21 -2.48 -1.92 soras . 2688 28951 26462 
22 =2.20 =iaee sare mers he a ~5926 25498 
23 =2.04 =i. Sa 2.0814 ea oo ~4326 65907 
24 -1.81 sk ae .53818 * £612 2950 26493 
25 i EL -1.48 «orbs 2606 93981 2.bo8e 
26 -1.42 =—1.47 {at ao .2644 ~4912 6433 
27 Ca al Ge =) ..493 5748 22526 894" 26476 
25 -1.62 -1.46 » 2679 -2661 . 1984 .6438 
29 sry a -1.47 - .asrs nave slaw 226° 
38 -.64 -1.46 » 4794 2693 » 1354 Pehl 
31 -.44 -1.47 41498 -2553 .2698 5 ee 
32 —,20 -1.47 24295 2658 2442 ~ag94 
33 -.85 -3. 98 * 5011 27468 » 1616 .6846 
34 21S “2.32 23562 22°18 . 1846 aor 
35 ~ 36 =1,¢%  < «Brees et Pig us73 ~5463 
36 .ao = ra | Pe Ye =] . 2636 ~8365 nals, 
ae ated =4 Ss - 5388 sevgl ke ar ~bS41 
33 o/s, -1.3¢ Pos he Jp he] s251s ~MS3 2.65954 
33 Awin, =4£.82 ~3902 wee 3s 28544 o blared 
48 i.cs -1.83 39983 = 29333 oy eeniteal 
+1 bao =1.5:3 ~2olG .25@1 W922 Part 3 =) 
42 Suse =) 55 »4458 25998 wa Aires oe) 
43 Say be -1.34 ~4421 Pa ee Sesitenesal 
+4 sao ee ON ee nee .2668 ~H99 weds 
45 cE bs. Sa |p A ae 4 ey 5 Slee ere 6481 


SL 


é 
s 
” 
- 
< 
_ 
. = 
rye 
7 * 
7 
a 
7 
7 - he 
1 a 
¢ 
a & 
° 
C= 





@ 
eo + 
e 
e+ 
c + 
= 
; 
= 
> 
é 
S 
a. 
* 
ait? 
ce 
; 
bé 


oe 


oe he 


‘ 


+ 
TeTirTy TTT 






ez és ih 


6 
Aa sees) 


beu4 








e 
oe 
7 
ewyvi i oO 


® 
2 © «¢ 


Ve a & 

‘ay 
ee » % 
= 


cs > 


z 
am 


TaA;T) 
ot me 
CS « 


2k 





TABLE VI. PROEE DATA, LOWER PLANE AT MIDSPAN (i = 5.3°) 


DATA FROM FILE LCBREDS:T1i4 
BLADE TO BLADE TRAVERSE MIDSPAN 


LOWER PLANE 


Point Loctin: Beta Ae Gttear Psibar Protilpar Se atar 
SREEEEREETEE TERT ELE TESLESESSPRAFERESSSPRE BEEPS REEKEFHEERELSRERSE TE EEKE 
1 -4.48 -42.45 293533 =, bare 28932 2345? 

a -3.98 -42.42 ~ 9603 -. 1444 83967 ~34r4 
3 -3.58 -42.42 2567 =, 1420 8956 22454 
4 -2.568 -42.44 9614 =< ae - 1631 34598 
s -2.00 -32.43 29741 -. 1469 0371 23524 
5 -1.58 -42.433 o 3767 =, 1493 .8375 ~ 8532 
7 -1.08 -42.43 SF 22 -.1481 24335 -some 
8 -.54 -42.43 ~ 9743 -. 1564 8921 3548 
3 8.88 -42.43 ~ 9799 = Looe W925 25548 
1a .38 =§2.42 . 3349 -.1532 ~8ss4 23585 
11 1.88 -42.43 »9814 -. 1684 8558 + Sots 
le 1.50 -42.43 9671 -.1628 ~@rrs .8635 
13 2.86 -42.43 «378? =, 1615 . 83998 » 8382 
is 2.56 -92.43 ~9742 ~, 1632 8971 23558 
15 3.04 ~42.41 me rae | -.1656 . 1859 . 3526 
i6 3.58 -42.42 . 3867 ~.,1663 8924 23589 
17 4.04 -42.44 ~9735 -,. 1631 .8968 . 33567 


32 
























1 : a a aé 

or 

a — ‘ aes 

i) 7 y r.€ i & -_—-—- 

in 

| : 

i 

i 

| | 

] i” a 'ot rb ee ie 

i °« e ‘ss* Sener egenege 

| - = 2¢ .ah= pid 

1 . , »* ae or.€=.. 7 

i . : : ve Be bout) 
: / v, b — oe © : 

i P ; : * 4 - é 1. 

. - f _- OF I~ - 

fi 





Sa3LE Vit. CENTER BLADE DATE: (4 = 5%3°) 


7C ca. Cpl Cpe Mach mine | 


FTE RLS RE RSS ACESS ECE HSER SHER ASR EEE e HERES SESH RRR E HE EEE ESHER EEE ESSE 


PRESSURE SIDE CENTER BLADE 


.! 


AGA? » 8054 26655 24573 21735 0°96 
»48108 -Ou19 26417 ~4148 1e49 3824 
»-9319 »48n0 23258 21999 ~2136 ae | 
.84035 -9112 23237 ~2882 2126 »O3946 
.9558 28215 24571 1327 2215 »8953 
21213 - 8383 .4818 ~ 1239 weacl 283995 
~1956 - 3452 -4768 ~leir 2e2oe? . 1004 
22695 -8576 2475? ~1128 ~2ed4 24999 
» 3433 » 8663 - +764 21133 »2e43 9998 
24192 ~87 15 .4892 -1366 22215 ~O958 
24938 -8736 «4871 e132? ~c2l3 20953 
-5659 -8727 ot7 71 1146 2241 8997 
25457 »007% 24956 1432 ~22u1 O97 
ie oT’ -9601 .4839 «1359 2c2ls »O0956 
»7384 .64e8r »4999 1469 2202 .9938 
- 8233 -O411 ~4+672 -8965 »e263 leer 
. 8633 6327 »3537 ~B°13 2232 Se 
«9632 -8238 24245 -9139 2354 ~ lad? 
239451 »-3123 ~3815 -. 85394 22943 18S6 

3550 » 9006 eerrs -.2482 ~2646 1175S 


SUCTION SITE CENTER BLADE 


-Olov 8227 =1. 3641 -3.2355 +97 o esd ra 
-9319 . 8316 -.f v6 72.1647 24249 « Teer 
23479 ~-4339 -. 4986 <1. oor, 2237S .ir Bs 
-3353 »45563 -. 3541 we Wee ee ie. »2675 16z2 
»1213 .8718 -.3431 =1.37 79 ~ 3639 »1615 
~ 1956 3978 ~_e re -1.2833 23535 1582 
22695 ~1179 -. 2598 -1.2684 odes » 1556 
23433 . 1309 —.203& , 71.1243 ~ 5453 21326 
4192 ~1399 -.1546 -1.9338 3377 1493 
»+9368 1432 - 28854 ~. 37098 2s2°8 2it4? 
~3659 1412 -. 86354 -.3599 ooeer ~i¢se 
2d? .1339 -4162 oF sd ~3ti6 ~ 1386 
»7 146 . 21269 »4631 - —. 6387 23425 1342 
2 334 ~ 1421 1456 -.4837 .2386 21298 
»3233 -4895 re i = tao 2 22385 21245 
3633 238755 -e3e6 -, 3363 2730 ~ Lete 
» 2432 »8593 -2636 -~.2748 2673 1135 
- 3431 24407 »2842 ~.¢560 2654 ee aS, 
3330 ~4206 »2948 gare ya - 2613 »ilet 
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TABLE VIII. ADJACENT BLADES DATA (i = 5.3°) 


See eg & Cpl Cpe Nach atte] 


| aa ie Rain te aE eet ee Se ie Mai ie Bee te te te Be eRe ee Rete Ree ee ee ee ee ee ee ee ee ee ed 


PRESSURE SIDE LEFT BLADE 


-12138 ~9383 »4014 -.8231 e402 1865 

~4132 28716 ~t519 ~9565 Zara 1612 

- 3233 »odil 244991 ~6E35 2302 1024 
SUCTION SIDE LEFT BLADE 

-1218 -67168 -.3438 -l.3r a8 3550 lels 

24192 - 1468 -. 1505 -1.6273 Say a 21431 

~3233 -3895 2197S —. 3943 2734 21240 
PRESSURE SIDE FIGHT BLADE 

»1218 » 8383 . 4648 . 8967 2278 .16148 

24192 -8716 »4558 28939 2256 . 1585 

. 3283 -O6411 24597 8329 2279 1814 
SUCTIGN SIDE RIGHT BLADE 

»1213 -87168 ~.3233-_ =1,.33513 3536 . Lede 

24192 - 14088 =, 1519 -1,82993 3374 1432 

.53233 »9395 -1911 -.4859 2385 1245 
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TABLE IX. BLADE PERFORMANCE DATA 


Present Results From Ref. 7 
8, 42.42 42.42 
Zz 5.3 5. 3 
85 1.85 0.4 
§ 10.44 9.0 
D 0.455 0.46 


E| 
oO 
. 

oO 
Ww 
~~] 
oO 
e 

oO 
rs 
ke 


3 
W COS 85 


—_-- 0.020 0.023 
20 cos By 
® cos 85 2 
(x 10°) 1.09 1242 
20 
AVDR 1.015 1.065 
‘ 0.351 
CPsraTic Sn 2 
mi -1.380 
Cw 1.385 
-0. -0.566 
Com 0.669 
ha =1.476 
, 1.330 
Cc -0.572 -0.645 
yB 
@, (in. 4,0) 25 
X .14 12 
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Fig. 4. 
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Photograph of Instrumented Blade 
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Fig. 8. Instrumented Blade 
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Fig. 11. Probe survey Data at Midspan of Lower Plane 
End Wails at 35°; Points iI -to~50 
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Fig. 12. Probe Survey Data at Midspan of Lower Plane 
End Wakls at 35°, Points 51 to 100 
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Fig. 13. Probe Survey Data at Midspan of Upper Plane 
End Wallis at 35°>.Points 1. to 50 


(PoLENUM ~ B) 1/0 of 


67 





\ 


























<a 





é 
a 
a 
‘ ‘he 





to alta 


ste vey 
y ae 
= 
4 
(6 - 


2 


za 


= } 
Paso 

iJ 

ad 


a 
, 

Q 

ur 











Piplenum) -P1 “Bred 
POINTS S1 TO 188 UPPER PLANE 35 DEG 


Fig. 14. Probe Survey Data at Midspan of Upper Plane 
End Walls at 35°, Points 51 to 100 
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Fig. 15. Probe Survey Data at. Midspan of Lower Plane 
End Walls at 30°, Points 1 to 50 
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Fig. 16. Probe Survey Data at Midspan of Lower Plane 





End Wallis at '30°., 
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Points 51 to 100 
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Fig. 17. Probe Survey Data at Midspan of Upper Plane 
End Walls at. 30°,.Points 1. to 50 
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Fig. 18. Probe Survey Data at Midspan of Upper Plane 
End Walls at 30°, Points 51 to 100 
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Fig. 19. Probe Survey Data at Midspan of Lower Plane 
End’ Walls at 50°, Points 1 to 50 
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Fig. 20. Probe Survey Data at Midspan of Lower Plane 
ma Walis at 50°> Peints* 51 to 100 
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Fig. 21. Probe Survey Data at Midspan of Upper Plane 
End Walls at 50° 
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Pig. 22. Probe Survey Data at Midspan at Lower Plane 
End Walls at 50°, Two Runs, Points 1 to 50 
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Fig. 23. Probe Survey Data at Midspan at Lower Plane 
End Wallis “at 50°, Ywo ‘Runs, Points Si‘to' 100 


(PSrENUM ~ Be I ede 


Toy 








Fen 14 be 






~ 
~ 
7 La 
~ . eo 
_ So ; 
7 
2 o : ag 
«3 > 
’ 4 
Pr} oo al 
of 7 - 
« . | * 
. Oa ey 
be eo“ ; 
“= 8% 7 
ry 7 - 4 » 
- ae + on 
=i : 












































Piplenum)-Pl-OQreaf 
* FIRST RUN 
+ SECOND RUN 

r 


= 
el rh F Ny ra 2, 
rN ‘, Pee er nf Rw =i wt 4 
+ \ + aie el mf ‘ed a a 7s Fr, : y ae 
Me " a) ¥ ‘oo: ¥ 
f 5 * Ay % VS 
[ 

4 eee ee eee 
os - be i ua ™ _ 
ats = | | | | | 

! 


INCHES 


Fig. 24. Probe Survey Data at Midspan at Lower Plane 
Ena Wails at 30°, Two Runs, Points 1° to 50 
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Fig. 25. Repetitive Samples with Fixed Probe Position 
(10" Left of CTR Midspan, End Walls 30°, 


Lower Plane (POL ENUM - Pamp) /Sree) 
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Fig. 26. Repetitive Samples with Fixed Probe Position 
(10" Left of CTR Midspan, End Walls at sore. 
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Fig. 27. Repetitive Samples with Fixed Probe Position 
(190” Left of CTR-Midspan; End Wallis at 30°, 
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Fig. 28. Repetitive Samples with Fixed Probe Position 
(on Centerline at Midspan, End Walls at 30°, 
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Fig. 29. Revetitive Samples with Fixed Probe Position 
(on Centerline at Midspan, End Walls at 3:0.%, 
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Pig. 33. Repetitive Samples with Fixed Probe Position 
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Fig. 34. Probe Survey Data at Midspan, Lower Plane 
16 Mesh Screen, Walls at 35°, Points 1 to 50 
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Fig. 35. Probe Survey Data at Midspan, Lower Plane 
16 Mesh Screen, Walls at 35°, Points 51 to 100 
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Fig. 36. Probe Survey Data at Midspan, Upper Plane 
16 Mesh Séereen, Walls at 35°, Points 1 \to, 50 


(P P,)/Q 


PLENUM ref 


90 











a a iby Fe 


a 
_ a 

tod) P= fed 

2. Behl Ay BE OT. BF 

: 7 


‘f acnd yeraed 
os 2 let . neers ge 
pag h tyt 


Cir 


Piplenum)-Pi-OQret 
POINTS S51 TO 166 UPPER PLANE 35 DEG TEMP SCREEN 





| 
| 
r+ 
«-* Pong en : geht 
ict pe “ — = a “ > ite to hee 
ay Oo : 
| 
.4 
| 
ee eee EEE | i ao 
re | cu ~- = = = 


[NtlHE'S 


Fig. 37. Probe Survey Data at Midspan, Upper Plane 
16 Mesh Screen, Walls at 35°, Points 51 to 100 
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Fig. 41. Probe Survey Data at Midspan, Lower Plane 
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Fig. 42. Probe Survey Data at Midspan, Lower Plane 
16) Meeh and 2, Mesh, Walls at 35°, Points 51:to 100 
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Fig. 43. Probe Survey Data at Midspan, Lower Plane 
4 Mesh Screen, Walls at 35°, Points 1 to 50 
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Fig. 44. Probe Survey Data at Midspan, Lower Plane 
4 Mesh Screen, Walls at 35°, Points 51 to 100 
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Pig. 45. Probe Survey Data at Midspan, Lower Plane 
5 Mesh Screen, Walls at 35°, Points 1 to 50 
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Fig. 46. Probe Survey Data at Midspan, Lower Plane 
5 Mesh Screen, Walls at 35°, Points 51 to 100 
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Fig. 48. Probe Survey Data at Upstream Midspan 
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Fig. 51. Probe Survey Data at Midspan 
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Probe Survey Data at Midspan 
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Fig. 54. Spanwise Probe Data Surveyed 1 in. from 
Suction Side of Centermost Blade 
ii = 5.3, (PorenoM P,)/Q), Upper Plane) 
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Fig. 55. Spanwise Probe Data Surveyed 1 in. from 
Suction Side of Centermost Blade 
(i = 5.3, X/X, Upper Plane) 
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Fig. 57. Spanwise Probe Data Surveyed 1 in. from 
Pressure Side of Centermost Blade 
(i= 5.3, .%/%, Upper, Plane) 
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Fig. 58. Resultant Blade Force Vectors by Momentum 
Balance (——--—--—--— ) and from Surface Pressure 
Integration (————) i= 5.3 
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Fig. 59. Measured Blade Surface Pressure Distribution 
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APPENDIX A 


MODIFICATION TO THE INLET GUIDE VANE SECTION 
OF THE SUBSONIC CASCADE WIND TUNNEL 


As discussed in Section I, Cina discovered during his 
test program that while the inlet flow to the test section 
was uniform in direction and uniform in wall static pres- 
sure, it contained a variation in velocity and stagnation 
pressure resulting from the wakes of the IGV's. Because the 
inlet guide vanes were spaced at intervals of two inches and 
the test cascade blades spaced at three inches, departures 
from strictly periodic conditions were detected from one 
test blade passage to another. 

To alleviate this problem, the inlet guide vane arrange- 
ment was modified so that the guide vanes were placed at l 
inch intervals. In order to preserve the option of reverting 
to a two inch IGV arrangement and because it was not possi- 
ble to machine the south wall to hold additional blades, a 
separate structure was placed between the bell mouth contrac- 
tion and the walls of the cascade. By mounting the IGV's 
in a separate unit which remained fixed once installed, 
hardware adjustments between tests associated with a change 
of end wall angle were greatly simplified. 

The new inlet guide vane assembly was constructed using 


two lengths of 10 inch steel channel as shown in Fig. A.l. 
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One set of guide vanes was mounted on the south side of the 
unit at 2 inch spacings. A second set of vanes was mounted 
at 2 inch spacings on the north side of the unit. When the 
unit was assembled the guide vanes mounted from alternate 
sides meshed, resulting in an inlet guide vane spacing of 

1 inch. The unit was provided with a single hand crank at 
the east end so that the vanes would be adjusted in unison. 
Once installed the complete structure could be left in place 
when the cascade north wall was removed to adjust air inlet 
angle. The one inch vane spacing ensured that periodicity 
at the test section would result for any test blade spacing 
which was a multiple of 1 inch. Equally important, the 
wakes remaining at the inlet to the test cascade would be 
greatly reduced as a result of closer spacing. 

Figure A.1l shows the details of the inlet guide vane 
unit, while Fig. A.2 shows the assembly in relation to the 
bellmouth contraction and the side walls. Figure A.3 shows 
the mechanism to adjust the inlet guide vanes. Figure A.4 
shows a view of the IGV assembly from the north side. A 
view of the Cascade Wind Tunnel partially assembled (north 


Side wall off) is shown in Fig. A.5. 
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Fig. A.i1. \Inlet Guide Vane 
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View of the IGV Adjustment Mechanism 
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Side View of the IGV Assembly 
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View of the Subsonic Cascade Wind Tunnel (North Wall Removed) 
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Fig. 






































APPENDIX B 


SELECTION AND INSTALLATION OF SCREEN MATERIAL 


Pankhurst and Holder [Ref. 17] show that the turbulence 
of an airstream can be increased by placing a coarse mesh 
across the flow upstream of the test section. One of the 
most effective methods which is used for the reduction of 
turbulence and non-uniformities also consists of placing a 
mesh screen across the tunnel. Screens used for this pur- 
pose are of a finer mesh and are placed at a greater dis- 
tance from the test section, and normally in the low-speed 
region upstream of the contraction in a conventional sub- 
sonic wind tunnel. By using such a screen the large scale 
eddies are removed at the expense of the introduction of a 
greater number of smaller eddies which decay rapidly. 

McEligot [Ref. 16] investigated the possibilities of 
reducing non-uniformities in the test section of the sub- 
sonic cascade wind tunnel. His investigation and recom- 
mendations were completed while the inlet guide vanes were 
‘Still at 2 inch spacings. McEligot concluded that some 
modification was necessary to achieve one percent uniformity 
for the mean velocity at the test cascade inlet plane and 
suggested several options. One of the options suggested was 
placing the turning vanes (inlet guide vanes) at a closer 


pitch. As explained in Appendix A, the pitch of the inlet 
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guide vanes was reduced from 2 inches to 1 inch. This new 
inlet guide vane arrangement did result in a one percent 
uniformity for the mean velocity at the test cascade inlet 
plane. 

The other approach suggested was the use of wire gauze 
screens. McEligot showed that the velocity distribution ap- 
peared to be largely dependent on a pressure drop coefficient 


K, defined by the equation 


Py = Pa 


kK = 
pv" 


where Py and P5 are the pressures upstream and downstream of 
the screen respectively. This pressure drop coefficient de- 
pends mainly on the blockage coefficient 8 defined by the 


equation 
8B = (1 - as)? 


where d is the diameter of wire used in the screen and 2% is 
the distance between the wires. This blockage coefficient 
is commonly referred to as "percent open area" in catalogues 
of industrial wire cloth and woven wire screens. 

For the velocities and flow angles used in the cascade 
wind tunnel, McEligot recommended using a wire gauze screen 
with a resistance coefficient, K, of 2.2 and a blockage 
coefficient, 8, of 0.47. However, since the new inlet guide 


vane arrangement resulted in a one percent uniformity for 
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the mean velocity and the pressure drop across a screen with 
a blockage coefficient of 0.47 was expected to be higher 
than could be tolerated for the desired test conditions, 
screens with a slightly higher blockage coefficient (higher 
percent of open area) were selected to be tested. 


The screens tested were of the following configurations: 


MESH WIRE DIAMETER (inches) BLOCKAGE COEFFICIENT 
~ .0410 -6989 
5 .0410 «6320 
16 .0105 ~s922 


Until the effectiveness of wire gauze screen in reducing 
non-uniformities in this facility was proven, a temporary 
installation of the screen material was considered adequate 
for testing purposes. The test screen was installed in the 
cascade wind tunnel by placing it between the inlet guide 
vane assembly and the north and south end walls. This ar- 
rangement placed the screen 7.25 inches downstream of the 
inlet guide vanes and19.3 inches upstream of the lower test 
plane. Figure B.1l shows the installation of the wire gauze 


screen. 
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APPENDIX C 


CASCADE PERFORMANCE PARAMETERS 


(by F. S. Cina; reproduced with minor changes from Ref. 7) 


The performance of a cascade is specified in terms of 
the deviation angle (4) and the loss coefficient (W) for 
given inlet conditions. In Ref. 1 the loss coefficient is 
shown to correlate in terms of the Diffusion Factor (D). In 
the present work, the performance parameters were calculated 
using the following expressions: 


1. Loss Coefficient (w) 


o = ——— (1) 


Pp P 
ty 1 


where the mass averaged pressure coefficients in Eq. (1) are 


O| 
\ 


defined in Appendix C of Ref. 7. It is shown in Appendix 
C of Ref. 7 that the effect of time dependent supply con- 
ditions are removed and the effect of AVDR is included in 
the use of Eq. (1). 


2. Diffusion Factor (D) 


WwW AW 
ms ie ae. 
a w. * 2ow (2) 


i it 
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3. Pressure Rise 


“eS (3) 
static Qa 


4. Blade Surface Pressure Coefficients 


os (4) 


a (5) 


5. Dimensionless Velocity 


- where V is the local velocity, ¥pe y2C Te is; the Tljmating” 


velocity and T, is the stagnation temperature. 
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APPENDIX D 


INSTRUCTIONS FOR PREPARING INPUT AND OPERATING 
QSONIC USING A RECTILINEAR CASCADE CONFIGURATION 


D.1 BACKGROUND INFORMATION 

QSONIC has the capability to calculate an axial, mixed 
or radial flow field and the test cascade can be rotating or 
stationary. The geometry of the streamsurface can be a 2D 
planar cascade or axisymmetric with varying channel thick- 
ness and radial position. The capabilities of QSONIC, be- 
yond those of previous cascade analysis methods (such as 
described in Ref. 8) include the ability to calculate 
through weak shocks with a peak relative Mach number less 
than 1.4, and completely around both leading and trailing 
edge regions of a blade profile. The blade shapes in the 
leading and trailing edge regions are not restricted to 
circular arcs. Detailed instructions for preparing input 
for a configuration other than an axial flow, stationary and 
rectilinear cascade may be found in Ref. 9. What follows 
are instructions for preparing the input applicable to the 
Rectilinear Cascade facility and running QSONIC on the Naval 
Postgraduate School's computer and associated operating sys- 
tem. It is assumed that the reader has a working knowledge 
of the NPS computer operating system and is familiar with 


Refs. 13 and 14. 
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QSONIC operates in two parts. The first part gener- 
ates a body (blade) centered mesh (geometry generation). 

The second part actually solves for the flow conditions at 
points in the mesh (flow solution). The data necessary to 
generate a mesh consists of a two-dimensional description of 
the blade shape. This is in the form of pairs of (X,Y) 
points on the surface, together with parameters that de- 
scribe the cascade layout, such as chord and stagger angle. 
Additionally, parameters describing the density of mesh 
lines complete the input for the geometry generation. 

For flow field calculations, the upstream flow condi- 
tions, convergence criteria and a schedule of meshes to be 
used should be input. If quasi-three dimensional effects 
are to be considered, a data file containing a description 
of the streamchannel's radial thickness and position as a 
function of distance along the stream surface is needed. 

For the case described herein, this was input by assuming a 
linear reduction in streamchannel thickness using a factor 
of 1/AVDR. This gave excellent results. (The output of 
another NASA code, Meridl [Ref. 15], can be used to input 
data to QSONIC for compressor flow field calculations. This 
program has recently become operational on the NPS computer.) 

The output of QSONIC consists of listings which con- 
tain an echo print of the input data, generated mesh coordi- 
nates on the blade surface, progress reports on the flow 


convergence and a list of the final velocities, pressures 
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and densities on the blade surface for each grid that was 


included in the schedule of solution meshes. 


D.2 INPUT DESCRIPTION 

The input for QSONIC falls into the following 
categories: 

Logical and case control parameters (NAMELIST PARAMS) 
Bulk data input: 

For geometry generation runs (NAMELIST INSTUFF) 

For flow solution runs: 

Mesh point storage files 

Streamchannel data file (for quasi-3D) 

The following is a description of the logical and case 
control parameters for the Rectilinear Cascade. Except for 
TITLE, the format for all these variables is in namelist 
form. The namelist is PARAMS. This information is taken 
from Ref. 9 and adapted to the Cascade Wind Tunnel. Name- 
list variables can be entered in any order. If a default 
value is listed, it is not necessary to enter that particu- 
lar variable. If the default value is listed as none, then 
a value for that parameter must be input. 

The following parameters apply to both the mesh gen- 
eration and the flow solution runs, but the values need not 


be the same. 
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Variable Default 


Name Type Value Description 
TITLE Alphameric None This is a one-line name for 


the case being run. TITLE 

Must appear on the first line 
of the data file that will be 
referred to as NAMELIST DATA. 


NOFLOW Logical False NOFLOW = .TRUE. if a run is 
to stop after generating a 
mesh, such as the mesh genera- 
f£ion° runs 


NOFLOW = .FALSE. for the flow 
solution run. 


MS Integer None MS is an array (max. dimension 
10) of values (max. value = 
25) for the number of grid 
lines in the mesh that will 
enclose the blade. 25 isa 
Satisfactory number for blades 
of solidity near unity. As 
solidity increases the maxi- 
mum value in MS should de- 
crease. 


NOZES Integer None NOZES is an array (max. 
dimension 10) of values (max. 
value = 49) for the number of 
grid lines in the mesh ra- 
diating from the blade on one 
surface. If a value of NOZES 
is greater than zero, a mesh 
with that many lines will be 
developed and stored in the 
file MESHGEN DATA. If the 
value is negative, QSONIC 
assumes that this MESH al- 
ready exists in the file 
MESHGEN DATA and will be read 
in. For Geometry Generation 
runs NOZES > 0 and for flow 
solution runs NOZES < 0. For 
electrostatic analog grid 
generator, NOZES must be odd. 
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BETA1 


BETA2 


GAMMA 


TOLS 


Real 


Real 


Real 


Real 


None 


None 


None 


BETA1 is the flow angle at the 
upstream boundary, in degrees. 
BETA] is measured from the aero- 
dynamic chordline to the direc- 
tion of flow; clockwise is 
negative. This can be obtained 
from the output of the cascade 
tunnel data reduction program 
"CX4431'. It is listed as "in- 
let air angle." Because of a 
difference in definition, it 

is necessary to use the usual 
inlet air angle minus the stag- 
ger angle for BETAI1. 


Example: The output of CX4431 
lists an inlet air angle of 
42.429°. The cascade is con- 
figured with a stagger angle of 
14227": 


42..429° 
- 14.27 


28.16 = BETAI1 


BETA2 is the flow angle at the 
downstream boundary in degrees. 
It 1s measured from the aero- 
dynamic chordline to direction 
of flow. Clockwise is negative. 
This can be obtained from the 
output of the cascade tunnel 
data reduction program 'CX4431'. 
It is listed as outlet air 
angle. Because of a difference 
in definition, it is necessary 
to use the outlet air angle 
minus the stagger angle for 
BETA2. 


This-.21s -bhe-catie.o0f specific 
heats. For the Subsonic Cas- 
cade Wind Tunnel the default 
value works well. 


This is an array of dimension 
10 of tolerances corresponding 
to MS and NOZES. Each grid 
solution will proceed until its 
TOLS value is satisfied. 
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MESH1 


MESHN 


LAMDAO 


CHORD 


Integer 


Integer 


Real 


Real 


Real 


None 


None 


None 


None 


This is an index in the 
arrays MS and NOZES of the 
first mesh to be generated 
and/or used for the flow so- 
lution. For geometry genera- 
tion runs, MESH1 selects one 
of the grids to be stored, 
provided NOZES(MESH1) > 0. 


This is an index in the 
arrays MS and NOZES of last 
grid to be calculated for 
geometry generation runs or 
used for the flow solution. 
For geometry generation runs, 
MESHN = MESH1. Subsequent 
flow solution runs then solve 
the case for all grids listed 
in MS and NOZES between index 
MESH1 and MESHN. For flow 
solution MESH1 normally 
points to the coarsest mesh. 


To insure equal spacing of 
grid lines over the entire 
mesh, 


NOZES(I) - 1 
KiI =— MESH 


Integer, 


where MESH1 < I S MESHN and 
senpe3 pit posy. Gath i2aF 
grid lines are doubled be- 
tween successive grids). 


Stagger angle of the blade 

row in degrees measured from 
the throughflow direction to 
the blade chord line. (Clock- 
wise is negative.) 


True chord of the blade, in 
the same units as the blade 
coordinates. 


Blade spacing in the same 


units as the blade 
coordinates. 
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The following parameters are required only for the 
geometry generation run. QSONIC gives the user a choice of 
two grid (mesh) generators. The Electrostatic Analog grid 
generator is applicable to any blade shape and any value of 
stagger or turning. The interpolation scheme grid genera- 
tor works for most blades except those with sharp leading 
edges where the edge radius is less than .5% of chord. The 
interpolation scheme allows the user to concentrate grid 
lines in areas of high interest around the blade. With the 
Electrostatic Analog grid generator no concentration of grid 
lines is available. For the cascade configuration used in 
this study, the interpolation scheme provided gross errors 
in the flow solutions, so the Electrostatic Analog grid 


generator was used with good results. 


Name Type Default Description 
NED Integer None Total number of body definition 


coordinates that are input. These 
are the X,Y points that describe 
the profile of the blade, with the 
first point repeated as the last 
point. 


KN Integer None KN is used to indicate which grid 
generator is to be used. KN = 0 
will call the electrostatic analog 
grid generator. For interpolation 
scheme, KN = number of body points 
on upper surface from the minimum 
to maximum X points, inclusive. 


NED and KN are required for either grid generator. If 


the electrostatic analog grid generator is used, no other 
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Parameters are used. 


If this grid generator is used the 


value from NOZES that is used must be odd. 


The following parameters are used only if the Interpo- 


lation Scheme is used to generate the mesh. 


Name 


RLE 


THETL 


THETT 


CAMPER 


STABAC 


Type 
Real 


Real 


Real 


Real 


Integer 


Real 


Default 


None 


None 


0.999 


Description 


RLE is the leading edge radius 
of the blade, with units the 
same as the coordinates de- 
fining the blade profile. 


RTE is the trailing edge radius 
of blade. 


THETL is the camber angle of 
the leading edge in degrees. 

It is measured from the aero- 
dynamic chord to the line tan- 
gent to the mean camber line at 
the leading edge; clockwise is 
positive. 


THETT is the camber angle at 
the trailing edge in degrees. 
It is measured from the chord 
line to the mean camber line at 
the trailing edge. Clockwise 
is positive. 


For blades whose chordline lies 
outside the blade profile, such 
as the DCA blading discussed in 
this report, extra grid lines 
Surrounding the blade are 
needed to interpolate the blade 
position. The truncated value 
of MS( )/CAMPER is added to 

MS( ). The maximum allowed 
value of MS( ) + MS( )/CAMPER 
is 30 grid lines. These are 
the grid lines that enclose the 
blade profile. 


STABAC is used only for test 
cases where no blade shape is 
to be input. Default value is 
usually adequate. 
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CHOP Real @,99 If leading edge or trailing 
edge radius is less than 2% of 
chord consult Ref. 9 for 
Clarification. Normally, 

RS: = CeCe St050.. 


SMOOTH Logical False For automatic addition of more 
blade definition points in the 
region of the leading and trail- 
ing edges set SMOOTH = .TRUE. 
This should be done for all 
blades except cusps and wedges. 


LEONLY Logical False LEONLY = .TRUE. for smoothing 
about the leading edge only; 
this is used if the trailing 
edge is a cusp or a wedge. 
SMOOTH must also be .TRUE. 


SLP1 Real 1 These parameters control the 
SLP2 Real 2 concentration of grid lines, if 
SLP3 Real 1 desired. The default values 
SLP4 Real 1 worked well for the DCA blades 


reported herein. For con- 
trolling the amount and loca- 
tion of the concentration 
consult Ref. 9. 


The following logical and case control parameters are 


required only for flow solution runs. 


Name Type Default Description 





MINF Real None This is the Mach number at the 
upstream boundary. This can be 
determined from the non- 
dimensional velocity X output 
from the cascade tunnel data 
reduction program 'Redd 5' and 
the relationship, 


ey x? 
y) M = eer: which yields 
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TCLS (I) is the tolerance for 
MS(I) and NOZES(I). There are 
three forms of input permitted. 


A) <-1.0 < TOES(Z) < 0.0: Cal- 
culations of the flow solu- 
tion will proceed until the 
relative circulation error 


N 
C catc ~ Cexacr 


Cexacr 


< |TOLS (I) | 


TOLS values between -1072 
and =~ DOV? are typical 
values for grids. This 
method of input is appro- 
priate only for lifting 
(non-symmetric blades) 
cases. 


By ond S0¢< ATOLS (IT) << b.0s - Cai= 
culations of the flow solu- 
tion will proceed until the 
average relative change in 
potential is less than the 
absolute value of TOLS(I). 


(2) < |TOLS (I) | 
> /avE 


Typical values should be 


between 10-”° and 107°. 

Cc). 1.0 .< SOL5(1).:.. Caloula- 
tions proceed until the 
number of iterations equals 
TOLS (I). 


Regardless of TOLS input, the 
solution for each grid will 
stop after 300 iterations if 
the TOLS criteria has not yet 
been made. All three forms 
were used for the cascade con- 
figuration reported herein 
with no discernible differences 
in results. 
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. FALSE. 


» FALSE. 


The default values of these 
parameters are adequate for 
flow conditions in the subsonic 
cascade wind tunnel. 


Number of iterations between 
intermediate printouts of re- 
Siduals and Mach number. The 
information controlled by this 
parameter is of limited value 
in comparing with measured 
data, so a value greater than 
10 reduces the amount of com- 
puter printout. For the study 
reported herein 40 was used. 


To list the flow quantities at 
all grid points in the last 
mesh set ALLOUT = .TRUE.. Un- 
less a very coarse grid is 
used, the output resulting from 
ALLOUT = .TRUE. would be ex- 
tremely voluminous and of 
limited value. Until the cas- 
cade is configured so it is 
possible to take data from 
between the blades, ALLOUT 
should be .FALSE.. 


QUASI3 = .TRUE. to activate 
streamchannel thickness and/or 
radius variations. The cascade 
wind tunnel has no radius 
variations, but to simulate 3-D 
effects the streamchannel 
thickness is reduced at the 
exit boundary by a factor of 
1/AVDR. This data is placed in 
a file used by QSONIC if a 
quasi 3-D solution is desired. 


This is the position of desired 
streamsurface data on the 
streamchannel file used if 
QUASI3 = .TRUE.. Currently the 
default value of 1 identifies 
the proper streamsurface data 
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in the streamchannel file. If 
the output from the NASA code 
"Meridl' is used for the 
streamchannel data, then by 
using different values of 
NSTRM, different streamsurface 
data may be used. 


RINF Real 4.0 This is the spanwise radius 
at the upstream boundary di- 
vided by aerodynamic chord. 
Radius effects are activated 
if RINF # 1.0. The current 
version of QSONIC allows the 
following cases. 


QUASI RINF Results 
{jy * FALSE: 1.0 Planar 2D flow 


2) 6. Tee. #1.0 Thickness on 
file; radius 
on file. 


Ss) a. let. 1.0 ‘Thickness on 
file; con- 
stant radius. 


Only tj) amd 3) “apply to >the 
cascade wind tunnel. 


WAKE Real 0.0 These parameters apply only if 
MINF2 Real 10.0 the test cascade is rotating 
OMEGA Real 0.0 and/or the downstream Mach is 
VAXIAL Real 999.0 Hear 1 .U. 

FLOCO Real 999.0 


At this point all of the logical and case control para- 
meters necessary to use QSONIC for the flow conditions pos- 
sible in the subsonic cascade wind tunnel have been discussed. 
The following is a description of the Bulk Data input for 


the geometry generation run and the flow solution run. 
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The format for all variables in the bulk data for mesh 
generation (geometry) is namelist form. The namelist is 


INSTUFF. 


Name Type Default Description 





H2 Complex None This is a table of points de- 
fining the blade profile. The 
real vart = X, and imaginary 
part = Y coordinate. The table 
begins at the point of maximum X 
value at the trailing edge and 
proceeds clockwise back around 
to the first point, which is re- 
peatec. The blade must be at the 
stagger angle and the origin at 
the point of minimum X. For 
Electrostatic Analog grids, the 
stagger angle must be positive 
(leading edge low, trailine edae 
high). The maximum number of 
points in H2 is 99 for the inter- 
polation scheme or 63 for the 
Electrostatic Analog. 


BUG2 Logical .FALSE. BUG2 = .TRUE. for a more detailed 
output of geometry generation. 
This will include the X,Y coordi- 
nates that define the mesh as 
well as second derivatives at 
grid points on the body. Except 
for trouble shooting this data is 
of limited value at the present 
time since there is no graphic 
output. 


The bulk data for flow solutions consists of a mesh 
file and the streamchannel data file if quasi-3D effects are 
to be calculated. The mesh file is created by QSOMNIC during 


the mesh (grid) generation run. No further inputs are re- 


quired from the user for the mesh file. 
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The streamchannel data file must contain a table of 
streamtube thicknesses, radial positions and corresponding 


X values along the streamsurface. 


Name Type Default Description 








CHO Real None CHO is the aerodynamic chord 
multiplied by the cosine of 
LAMDAO. (LAMDAO = stagger angle) 


NRSP Integer None NRSP is the total number of data 
points in each of the tables of 
thickness, radial position, and X 
location. If MRSP is 2, a linear 
distribution is obtained between 
the endpoints given. NRSP = 2 
was used for the study reported 
herein with good results. 


RM Real None Array of corresponding X loca- 
tions for thickness and radius 
data values. xX = 0 represents 
the leading edge of blade, with 
the blade at stagger angle. The 
units can be any consistent 
length scale common to CHO, RM, 
RMSP and BESP. Inches were used 
in this study. 


RMSP Real None Spanwise radial positions of 
streamsurface at the X locations 
given in RM. RMSP was not used 
in the current study. 


BESP Real None Array of streamtube thickness 
values at the X locations speci- 
fied in RM. For the study re- 
ported herein, at X=Oa 
streamtube thickness of 1.0 was 
arbitrarily selected. The 
streamtube contraction through 
the test section was similated 
by reducing the thickness at 
X= 0, by a tactor of L/AVDR at 
the trailing edge. Duval [Ref. 
3] explains AVDR. 
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D.3 PREPARING INPUT FILES 


QSONIC was originally configured to use several input/ 


output devices while reading data and generating output. 


The input/output devices are listed below as used by QSONIC. 


I/O Unit 


2 


Eb 


18 


23 


Usace 


File containing streamchannel data. This is used 
only if QUASI = .TRUE.. 


Standard card input; NAMELISTS PARAMS, INSTUF. 
Standard printed output. 


For mesh generation runs, coordinates of all mesh 
points are written here. For flow solutions, X, 
Y, velocities, pressures and MACE are recorded 

for graphic display. The program currently has no 
graphic output capability. No user action is 
necessary to create this file. 


Used as temporary storage. No output is stored 
here. No user action is necessary in conjunction 
with this file. 


Previously developed mesh coordinates are read in 
from this file during the flow solution. After a 
mesh generation run, the user must create this 
file and put in it the data from I/O unit 13, so 
that during the flow solution QSONIC can read in 
the mesh points. 


OSONIC is presently configured to operate with the CMS 


system of the IBM 370 computer. This system provides a high 


degree of flexibility in parameter selection. With this 


system, all the input/output units previously mentioned are 


on the disk space assigned to the Turbopropulsion Laboratory 


(TPL). 


Access to QSONIC and the TPL disk space can be ob- 


tained through the Director of TPL. 
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The first step in using QSONIC is the creation of the 
data file necessary for the mesh generation run. This is 
done using the standard procedures of the XEDIT function of 
the CMS operating system. Reference 14 has specific in- 
structions for creating new files. The filename and file- 
type for the data used in this study was NAMELIST GFOMD. 
Once the data file is opened, the necessary data is input 
beginning at column 2 of the virtual card. Since the varia- 
bles are in namelist form, they can be input in any order. 

Table D.1l is an example of the data file necessary for 
the mesh generation run. The TITLE must appear on the first 
line (FORMAT = 20A4). After the TITLE, the logical and case 
control parameters are input after the namelist GPARAs. 
When all the case control parameters required are input the 
PARAMS namelist is closed with SEND. On the next line of 
the data file the bulk data for the mesh generation run is 
input in namelist form, with the namelist &INSTUF. The 
Hl = 100*(0.0,0.0) that appears after &INSTUF on Table D.1 
was used on earlier versions of QSONIC, but is not used in 
the present version. It should, however, appear in the data 
file before the H2 variables (X,Y coordinates defining the 
blade profile). 

At this point, some Giscussion of the coordinates de- 
fining the blade profile is warranted. Table D.2 is a 
listing of the X and Y coordinates of the DCA blading used 


in this study. Figure D.1 is a plot of these coordinates. 
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Recall that the coordinates defining the blade profile for 
QSONIC must be for the blade at the stagger angle and mini- 
mum X at the origin. The coordinates of Table D.2 were 
translated and rotated using a coordinate transformation 
routine for the HP-67 programmable calculator. These new 
coordinates appear in the namelist INSTUFF on Table D.1. 
Figure D.2 is a plot of the translated and rotated coordi- 
nates. It is highly recommended that such a plot be made 
for any new blade profiles, to ensure that the original 
coordinates are translated and rotated properly. 

The second step in using QSCNIC is the creation of the 
Gata files necessary for the flow solution run. The file 
used for the flow solution in this report is on the ™PL disk 
space with a filename/filetype of NAMELIST FLOWD. The sim- 
plest way to open this data file is to use the YEDIT func- 
tion, as discussed a Ref. 14, to start a new file. Then 
input the same data as is in the data file for the mesh 
generation run usino the X=DIT subcommand GF™ (filenane) 
(filetype). The appropriate changes anc additions can then 
be made to this file. Table D.3 is an example of the data 
file just discussed. 

Two more data files are required for the flow solution. 
The data for one of these is created by the mesh generation 
run. The other file contains the streamchannel thickness 


data for implementing quasi-3D effects. 
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After the mesh generation run, a file with the file- 
name/filetype MESHGEN DATA will appear on the disk. Create 
a new file with the filename/filetyne MESHIN DATA. This is 
most easily done by issuing the command 'XEDIT MESHIN DATA'; 
then use ‘GET MESHGEN DATA'. This file contains the pre- 
viously developed mesh coordinates. 

The streamchannel data file should have the filename/ 
filetype of DATA5D DATA. The format for the data file is 


shown below. 


Virtual Card Column No. Variable Name 
Z BLANK 
2 BLANK 
3 21-30 CHO 
: BLANK 
5 36-40 NRSP 
6 BLANK 
7 BLANK 
8 
9 
10 
11 
12 1-80 (8F10.5) RM 
As veel 
1 1-80 (8F10.5) RMSP (not used in 
this study) 
As needed 
1-80 (8F10.5) BESP 


Table D.4 is an example of the data file used in the 


present study. Since NRSP = 2 was used for this study a 
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linear distribution is assumed for the streamtube thickness 
values and only 2 values of RM and BESP are required; there- 


fore, only 1 virtual card was required for each array. 


D.4 PROGRAM OUTPUT 

The output generated by QSONIC for the geometry genera- 
tion run includes a printed listing (I/O unit 6) and a mesh 
point file (I/O unit 13). The printed listing under the CMS 
system I/O unit 6 is normally the computer terminal unless 
the command 'FILEDEF 06 PRINTER' has been invoked. It is 
unusual for the program to run properly the first time, so 
initially it is helpful to have the printed listing appear 
at the terminal. Once the program is running properly the 
output should be sent to the line printer. 

The flow solution run output consists of a printed 
listing (I/O unit 6) and a plot data save file (I/O unit 
13). Once the flow solution is running properly the printed 
listing should be sent to the line printer. 

Table D.5 is an example of the output generated by the 
geometry generation run. Figure D.3 is a plot of the grid 
output points on the blade surface, horizontal chord, pro- 
duced by the mesh generation run. Figure D.4 is a plot of 
the grid output points with the blade at the stagger angle. 
Table D.6 is an example of the output generated by the flow 


solution. 
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A detailed explanation of the printed output for the 


program QSONIC may be found in Ref. 9. 


D.5 RUNNING THE PROGRAM 
The files on the TPL disk space that apply to QSONIC 

are listed below: 

QSONIC EXEC Al 

QSONIC FORTRAN Al 

QSONIC TEXT Al 

NAMELIST GEOM Al 

NAMELIST FLOW Al 

NAMELIST GEOMD Al 

NAMELIST FLOWD Al 

DATA5 DATA Al 

DATA5D DATA Al 

QSONIC EXEC sets the input/output devices required to 

read and store data. QSONIC FORTRAN is the source program. 
To document the changes to QSONIC necessary to use the code 
with the IBM 370 operating system and serve as a reference 
for future users, a program listing is included at the end 
of this appendix. QSONIC TEXT is the computer executable 
code created when QSONIC FORTRAN is compiled. NAMELIST GEOM 
and NAMELIST FLOW are the data files for the geometry gen- 
eration and flow generation respectively for the example in 
Ref. 9. DATA5S DATA is the streamchannel data required for 


the quasi-3D solution for the example in Ref. 9. 
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NAMELIST GEOMD is the data file for the geometry gen- 
eration for the DCA blading used in the studv reported 
herein. NAMELIST FLOWD is the file for the flow solution 
for the cascade configuration used in this study. 

QSONIC expects the input data to be in a file on the 
TPL disk space named NAMELIST DATA. Since the first time 
QSONIC is run is to develop the body centered mesh, the file 
NAMELIST GEOMD must be renamed NAMELIST DATA, using proce- 
dures specified in Ref. 13. Because QSONIC requires large 
amounts of virtual memory, extra storage must be defined for 
the code to operate. This is accomplished by issuing the 
command ‘DEFINE STORAGE 1504K'. 

With the data file renamed and more storage defined, 
type 'QSONIC' to load the program. The output will appear 
on the terminal screen unless FILEDEF $6 PRINTER was invoked 
prior to loading the program. 

After the mesh generation is complete, rename NAMELIST 
DATA to NAMELIST GEOMD and change NAMELIST FLOWD to NAMELIST 
DATA. Create a data file with the filename/filetype MESHIN 
DATA. The elements of this file are the same as the elements 
in the file MESHGEN DATA that was created by the mesh gen- 
eration run. The necessary input/output files are now con- 
figured for a flow solution run. Issue the command 'OQSONIC' 
to begin execution. 

If the program output appears at the terminal it is 


possible to have some I/O error messages appear with the 
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output. This is because the write statements in QSONIC are 
formatted for the 132 character long line of the printer. 
These errors do not affect the validity of the program 
output. 

The explanation for any error or condition message 


generated by QSONIC can be found in Ref. 9. 


D.6 QSONIC UPDATE 

Recently an improved version of QSONIC was reported by 
NASA Lewis Research Center [Ref. 13]. The new version re- 
quires less virtual memory and executes approximately 303% 
faster than the version presently in use at NPS. Also, the 
output appears in a different format than is described in 
this appendix. Reference 18 describes the most recent ver- 


sion of QSONIC in detail. 
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INPUT DATA FOR MESH GENERATION RUN OF QSONIC 
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TEST BLADE COORDINATES 
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INPUT DATA FOR FLOW SOLUTION RUN OF QSONIC 


TABLE D.3. 


7 

e 

oO 

e 

N 

- 
~_ *- 
ooe# 
Ont 
faa) ous 
. orar 
Wwe Ow 
a a >>wW 
Tie oe 
ern ~O 
NW ea * 
TWO NO su 
- 3 ee a he | 
we wd 
CHIN thik 
Zz e000 * 

Oe oy mu 


WUIMPN ext oo 
reat i VOU O 


Lith *ia 
Is eVect © 
OQONlOwrstwWe 
BsZzWnNnace 

O eZ eZ il 
InrTocr iw 
OLRM x 
<0 swuL sna 
ZZOzZzt + = 

— oe se 1O 

WY esters « 

=wWi 

ad -+OacwoO ii 


ANNO Se S22 Ii i 


u + 


—< 

@o 

= uw 

Wo + 

FN oO 

cor 9 

Oo i=) 
Se Rd = 
oow oO 
a ou 
- - =O or 
=) ooe =f ° 
+ e#+O OO 2D 


NWW WW OP NW 9 + WW 
SOD MNMENAON eMWUIDORANNWHY ST 
MANA AMO © « eDOHODHDOOMNOY 
SHON TOSOD NM TORMNAMED 

NIN Set ww 8 HD INOT ETON 
Roc secinaagsceme ee 0Q HOM ee tet tetas 


ees ce eG eomm—m © -t @ 2 oe ees @ 
Socccocane *} © eOODCOCOOo 
+ NO eS WMO 


Wee ee erfOMOoMe ee & & & & & oe 
LO Stet Setet§ DELON NO ott ett ett HO 
BOOSCSO «# «¢ MWHiNmDOOOCODCOCO e« 
OFF He FOO 2 § HOH HHH HHH HO 
CDBWW WWW 9000 WW WWW WLU Ly 
PFN ORD & © mee wwe OE NNN SIE A o 
emt OPO Sete & & PDO ORSON MO 
CRN OMNNOCOCOC mem eRmOM~R OOM ¢ 
~NO VI + + FO RFOOMPNVOMNGAS DO 
RPT OMNWWWRON SING OMS 
—~ ee eo OPO lM & © 0 © © 8 oe oO 
AO COOP AMNNWMHOOCCOOCOOCOMN 
CO wee were OPO OM OO we wwe wr re ef 
> eee er ej QOOOO+ we we we mB ee ew 
US ee ee OY ed et N Uo eee oe 
Ott mtn 8 8 8 0 ON et ted ted ee tt 
SOC OO0 COO Cet OO0D0COCCOCOCO 
De ee eww OM $F H HH HHH H+ 
NWWWU ® & we of et WL LULU LU ls 
ASANOWIDD ET OMAN DANAE HO 
2OM STH ORMItD 0 8 8 ON PMNRODR 
OHA DNOMDAPO OOO OSIM SNS H1.0 
NN SSO Dt Ow wr ONMT ETT OO 
eset ested OO) & © & wetted edt et ted ted 
—~—t 8 8 8 8 8 8 e taomawwe—m 8 8 88 8 8 8 8 8 
COCO OCOCC0CFC0Os00rCOCCCCCOoO 
+ ~ OU 
OW eke we we wm we POOFO ke we me wm we ee 
oT et tet tet SHO ON Steet et ted 
ODNOOGTOOCOO «+ e&  eOCCOQOOCOO 
Nee ee HHH AO HOO +H HH HHH HH 
COwWWWWWtU O WLW WW WL LU 
wt OM AN DMN e ] oe EME MNNO Sh 


156 


PNAIDM HAD MOUND Wie BHP OOM MeN SLALL 
3 = DOOSDNDMOONMM O0OT OD HNOUVM aN 
OE NOMAD Nt NODWND SOMEAIm® Il 

SOT eS LON 4 LPO Net tO ORD SHINO OST TINO 
ose eot?eeoeeevreetr eee easeeese MDZ 
AONW I PD Wr HOO COO DOC OCCO OCC OCOOCCOCCCCOCOnNWw 

















— Mele A Py hare tas (OOO te, 


- 


Pr et » 3Qe .-.,@ Sooo 


i 
7 


t coer SE SSs Salome 

























ve VAG OUBLEL 


, aA Or ilies wot oma) “ays ne eae 
. ao te «© &ret Gave 
metiwy POoOOOhAAOe. 
Shes Bre 6 ewe ter 4 
2 a Ee (ee, Qe ote Fee 
_ ~ >) * es ih Oe mn 
TIO OO [ew > oO Ooc 
Orrng) G4 - 
*»** © &©@ © @ & VAD >+t"<<<6 iW 
in os * 1-4) wy @eeeene Bete fy 
a » mine « 6 eooaosS 
2 @ ¢ e+ OCSCed « © 
ie [Oc liad oo 
* - oan = > aye ge 
Mises 4 ik «© @ tees Cb TP Gee a 
prof yo» er DDS CD | 
a bee) Se 7 oS ee Oe 
; ad i Pde Mat re ® gas 
* 2 6 @ a Pe bese e+e 


-OoSe2 Gin ate A BID 00—% 

~~ awe + he \ } Caeerns amet c 

e «4 @ oo & co C$oOOvl..e-« «es 6 7 
ome," a ee 

iT Cmmexyim lithe » = © & © R49 beeps 

SODOOSS Gah “MO oO ORCRon 7 
**s ©4648 6 UO Prune & ab > eS & 

‘AC! MT ade be be oe ee LU AA 


SSM Ne © © 6 enh ep eteege 
a Narr a Tish lia open 
CO DO Cy ere ar BD a 
0b Py, ae ae s dlind FE COFD vlna tens mee aes 

* ‘tt! hina 6 4 6% 6 6. 
OG FOS ASO OOIOOO SO 
wri “ 


* - @ & @&@ arnt} On «=< 6 9 
meme eee ew @ ee, Ar Ping retin 













(24% 6 42COO4 446086 5 
er uantetAmD oS a UA a 7 
6 LEAR A Ie (em, hei A AG 


pd ta WER Oo 





oer 
~ = 


TABLE D.4. DATA FILE FOR QUASI-3D SOLUTION 


STREAMLINE DATA FOR DCAI1 
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SAMPLE OUTPUT FROM MESH 


TABLE D.5. 
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